%% FigureS20A positive side
% i = 22,32,42,52
for i = 32

T = i;

load(strcat('FigS20\Q00',num2str(T),'Icpos.mat'));

load(strcat('FigS20\Q00',num2str(T),'Icposh.mat'));

figure(233)
Ic1 = histogram(Ic5Kpos./0.0016,[300:3:1000]); hold on;
Ic2 = histogram(Ic5Kposh./0.0016,[300:3:1000]);


hold off

end

xlim([500 1000]);
ylim([0 150]);
hold off
%% FigureS20A negative side
% i = 22,32,42,52
for i = 22

T = i;



load(strcat('FigS20\Q00',num2str(T),'Icneg.mat'));

load(strcat('FigS20\Q00',num2str(T),'Icnegh.mat'));

figure(234)
Ic1 = histogram(Ic5Kneg./0.0016,[-1000:3:-300]); hold on;
Ic2 = histogram(Ic5Knegh./0.0016,[-1000:3:-300]);


hold off

end

xlim([-1000 -500]);
ylim([0,150]);
hold off

%% FigureS20B
figure(12)
T = 13 % T = 13,14,15,16 
load(strcat('FigS20\O00',num2str(T),'.mat'));
plot((1:40001)./2000,datatest(1:40001,1)./0.0016,'color',[0 0.4470 0.7410],'LineWidth',1);hold on;
plot((1:40001)./2000,datatest(1:40001,2)./(0.0016*4),'LineWidth',1);
hold off;
xlim([0 20]);ylim([-1400 1400])

%% FigureS20C
twst = [43.7,44.5,40,43.5,41,40.5,39.5,30,34.8];
deltapos = [7,4.5,7,2,0,1,59,4,11];
errpos = [5.2,5.2,13.8,5.9,5.8,3.2,23.5,4.9,14.0];

deltaneg = [0,4.1,62,6,29,15,46,6,22];
errneg = [5.2,4.8,14.3,6.4,23.8,3.3,12.3,5.2,21.7];

figure(942)
h1 = errorbar(twst,deltapos,errpos,'o','DisplayName','|J_{c+}^{Full}-J_{c+}^{Half}|');set(h1, 'MarkerFaceColor', get(h1,'Color'));hold on; h2 = errorbar(twst, deltaneg,errneg,'o','DisplayName','|J_{c-}^{Full}-J_{c-}^{Half}|');set(h2, 'MarkerFaceColor', get(h2,'Color'));hold off;
yticks([0:100:300]);
xticks([25:5:45]);
xlim([25 45]);
ylim([0,300])
legend

